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Context. We present a Chandra analysis of the X-ray spectra of 56 clusters of galaxies at z ^ 0.3, which cover a temperature range of 
O 3 < kT < 15 keV. 

\& Aims. Our analysis is aimed at measuring the iron abundance in the intra-cluster medium (ICM) out to the highest redshift probed to date. 

Methods. We made use of combined spectral analysis performed over five redshift bins at 0.3 ^ z $ 1.3 to estimate the average emission 
(-H weighted iron abundance. We applied non-parametric statistics to assess correlations between temperature, metallicity, and redshift. 
Q | Results. We find that the emission-weighted iron abundance measured within (0.15 - 0.3)R rjr in clusters below 5 keV is, on average, a factor 
of ~ 2 higher than in hotter clusters, following Z{T) ^ 0.88 T~ 0A1 Z G , which confirms the trend seen in local samples. We also find a constant 
average iron abundance Z Fe ^ 0.25 Z Q as a function of redshift, but only for clusters at z ~ 0.5. The emission-weighted iron abundance is 
significantly higher (Z Fe = 0.4 Z Q ) in the redshift range z — 0.3 - 0.5, approaching the value measured locally in the inner 0.15i? v ; r radii for a 
mix of cool-core and non cool-core clusters in the redshift range 0.1 < z < 0.3. The decrease in metallicity with redshift can be parametrized by 
a power law of the form ~ (1 + z)~'" 25 . We tested our results against selection effects and the possible evolution in the occurrence of metallicity 
and temperature gradients in our sample, and we do not find any evidence of a significant bias associated to these effects. 
Conclusions. The observed evolution implies that the average iron content of the ICM at the present epoch is a factor of ~ 2 larger than at 
z — 1.2. We confirm that the ICM is already significantly enriched (Z Fe = 0.25 Z G ) at a look-back time of 9 Gyr. Our data provide significant 
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d ' constraints on the time scales and physical processes that drive the chemical enrichment of the ICM. 
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1. Introduction 

Clusters of galaxies are the largest virialized structures in the 
Universe arising from the gravitational collaps e of rare high 
peaks of primordial density perturbations (e.g. Peebles 1993k 
Coles & Lucchinl l 1995HPeacockHl999r.lRosati et al]l2002t IVortl 



20051) . As a result of adiabatic compression and shocks gener- 



ated by supersonic motion during shell crossing and virializa- 
tion, a hot thin gas permeating the cluster's gravitational poten- 
tial well is formed. Typically this gas, which is enriched with 
metals ejected form Supernovae (SNe) e xplosions through sub- 
seque nt episodes of star formation (e.g. lMatteucci & Vettolani 
Il988l) . reaches temperatures of several 10 7 K and therefore 
emits mainly via thermal bremsstrahlung in the X-rays. At such 
temperatures most of the elements are either fully ionized or are 
in a high ionization state. Strong emission lines may originate 
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by collisional excitation of K- and L-shell transitions in highly 
ionized elements, such as H- and He-like iron, oxygen, silicon 
or sulfur. In the isothermal approximation, the line intensities 
depend on the abundances of heavy elements, while the contin- 
uum intensity is mainly due to hydrogen and helium. Therefore 
the equivalent width of a line, under the reasonable assumption 
of collisional equilibrium, gives a direct measurement of the 
abundance of the corresponding element. 

The chemical enrichment of the intra-cluster medium 
(ICM) is an unambiguous signature of star formation in cluster 
galaxies. A knowledge of the history of the ICM metal enrich- 
ment is also necessary in understanding the mode and epoch of 
cluster formation and the thermodynamic evolution of the clus- 
ter baryons in their hot and cold phases. In this respect, measur- 
ing the properties of the ICM at high redshift is important for 
constraining the physical processes involved in the diffusion of 
energy and metals within clusters. In particular, the chemical 
evolution of the ICM sets constraints on the SNe rate in the 
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cluster galaxies (see Ettori 2005). The SNe explosions are, in 
fact, the main contributor to the metal enrichment and are also 



expected to provide a sourc e of ICM heating (e.g. lPipino et al 
l2002tlTornatore et alJliooj) . 

In addition, much emphasis has recently been given to the 
evolution of global scaling relations for the ICM, such as the 
luminosit y-temperature relation and the entropy-temperature 
relation fa olden et ai1l2002l; Ivikhlinin et alJl2002t lEttori et al 



2004; Pratt et al. 2006; Maughan et al. 2006). Such scaling re 



lations represent a further signature of galaxy formation and 
super-massive black hole activity, whose relative contribution 
to th e ICM energetic i s under investigation. 

Tozzi et al.l (120031 hereafter Paper I) measured for the first 



time the average iron abundance in the ICM of hot clusters 
of galaxies out to redshift z — 1.3. Their findings, which 
mostly probed the redshift range 0.5 < z < 1.3 (with 90% 
of their sample in this range), suggested that the mean iron 
content of the ICM is approximately constant with a value 
Zf e a 0.3 zJ3 for clusters with temperatures kT > 5 keV. 
When comparing their results at high redshift with the local 
values of the iron abundance, they referred to the typi cal value 
Zp e ^ 0.3 Z available from the lit erature at that time ( Renzinil 
1997tlDe Grandi & Molendill200lh . Therefore, they concluded 



that no evolution of iron abundance with redshift was found out 
to the probed redshift. 

More recently, the determination of an average local value 
of the ICM iron content has become a harder task due to the 
different metal distribution inside cool-core and non cool-core 



clusters ( De Grandi et al. 20041) . The former show a central 
peak of iron abundance with Zf e 0.6 - 0.8 Z Q and a plateau 



at Zf> 0.3 Zq in th e outer regions (see iTamura et al.l 12004 



Vikhlinin et al. 2005), and the latter a somewhat lower value 



(Zf e - 0.2 - 0.3 Z Q ) constant with radi us. The iron-rich cores 
typically have a size of 100 kpc (see |Pe Grandi et al.l 12004; 



Vikh linin et alj 120051) . Therefore, it is important to take into 



account any effects due to different physical apertures when 
comparing different samples from the literature. 

Studies of local cluster samples have also found an in- 
crease in the ir on abundance in clusters with temperatures 
£ 5 keV (see Arnaud e t al.l Il992t iFinoguenov et alJ l2001t 



Baumgartner et al.1 



2005), whose physical interpretation is still 



a matter of debate. 

Single observations of high-z (z 



> 1) clusters confirmed 



thatZ 



2004 



0.3 Z B or higher is common in the ICM (Rosati et al 



Hashimoto et 



higher is 
alJ2004 



This implies that the last episode 
of star formation in clusters of galaxies must have taken place 
at earlier epochs in order to significantly enrich the diffuse 
medium with metals. Therefore, the study of the iron abun- 
dance at high redshift is expected to place strong constraints 
on the star formation history of cluster galaxies and on its ef- 
fects on the thermodynamics of the ICM. 

Here we present a significantly improved analysis com- 
pared to our previous work by substantially extending the sam- 



1 Solar abundan ce values are set to those provided by 

lAnders & Grevessel dl989l) ; in particular, the solar abundance of 
iron atoms relative to hydrogen is 4.68 x 10~ 5 . See Sect. 3.1 for 
discussion of the solar iron abundance. 



pie (56 clusters instead of the 19 presented in Paper I) and by 
using the most recent Chandra calibrations in order to have an 
up-to-date data reduction at the time of writing. This increase in 
statistics allows us to investigate the relation between the iron 
abundance and global temperature of the ICM at high redshift 
and to derive a more robust measurement of the cosmic evolu- 
tion of the average iron abundance in the ICM, which is then 
compared with predictions from the cosmic star formation rate. 

The plan of the paper is as follows. In Sect. 2 we describe 
the data reduction procedure. In Sect. 3 we describe our spec- 
tral analysis (Sect. 3.1) and present the main results, subdi- 
vided as follows: single source analysis (Sect. 3.2), correla- 
tion between iron abundance and temperature (Sect. 3.3), evo- 
lution of the average iron abundance as a function of redshift 
(Sect. 3.4), and proper comparison with the local iron abun- 
dance (Sect. 3.5). In Sect. 4 we discuss the implications of our 
findings and in Sect. 5 we summarize our conclusions. In the 
Appendix, we describe in detail the spectral simulations per- 
formed to investigate the possible spectral-fitting biases due to 
unresolved temperature and metallicity gradients in our sam- 
ple. We adopt a cosmological model with Ho = 70 km/s/Mpc, 
Q.M = 0.3, and £_a = 0.7 throughout. Quoted confidence inter- 
vals are lcr unless otherwise stated. 

2. Sample selection and data reduction 

2. 1. Chandra data 

In Table [T] we present the list of Chandra observations ana- 
lyzed in this paper. The selected sample consists of all the pub- 
lic Chandra archived observations of clusters with z > 0.4 as 
of June 2004, including 9 clusters with 0.3 < z < 0.4. Some 
of them were already presented in Paper I. Data reduction is 
performed using the CIAO 3.2 software package with a recent 
version of the Calibration Database (CALDB 3.0.0) including 
the correction for the degraded effective area of ACIS-I chips 
due to material accumulated on the ACIS optical blocking fil- 
ter at the epoch of the observation. We also applied the time- 
dependent gain correctiorQ, which is necessary to adjust the 
"effective gains", which have been drifting with time due to an 
increasing charge transfer inefficiency (CTI). Most of the ob- 
servations were carried out with the ACIS-I instrument, while 
for some clusters (see Table[T|) the Back Illuminated S3 chip of 
ACIS-S was also used. 

We started to process data from the level=l event file. 
For observations taken in the VFAINT mode, we run the tool 
acis.process.events to flag probable background events 
using all the information of the pulse heights in a 5 X 5 event 
island (as opposed to a 3x3 event island recorded in the FAINT 
mode) to help in distinguishing between genuine X-ray events 
and artificial events that are most likely associated with cosmic 
rays. With this procedure, the ACIS particle background can 
be significantly reduced compared to the standard grade selec- 
tion^. Real X-ray photons are hardly affected by such cleaning 
(only less than 2% of them are rejected, independent of the en- 
ergy band, provided there is no pileup). We also applied the 



2 http://asc.harvard.edu/ciao/threads/acistimegain/ 

3 http://asc.harvard.edu/cal/Links/Acis/acis/Cal_prods/vfbkgrnd/ 
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Cluster 


z 


Obs. Id." 


Exp. [ks] b 


Mode c 


R e * ["] d 


Net Counts c 


MS 1008.1 - 1224 


0.306 


926 


44 


I-V 


108.5 


9260 


MS 2137.3 - 2353 


0.313 


928 


33 


s-v 


79.0 


33000 


Abell 1995 


0.319 


906 


56.4 


S-F 


103.0 


30100 


MACS J0308.9 + 2645 


0.324 


3268 


24.4 


I-V 


123.0 


11200 


ZwCl 1358.1 +6245 


0.328 


516 


48.3 


S-F 


88.5 


19800 


MACS J0404.6+ 1109 


0.355 


3269 


21.6 


I-V 


157.0 


3100 


RXJ0027.6 + 2616 


0.367 


3249 


9.8 


I-V 


108.0 


960 


MACS J1720.2 + 3536 


0.391 


3280 


20.8 


I-V 


103.0 


6670 


ZwCl 0024.0+ 1652 


0.395 


929 


39.5 


S-F 


64.0 


3150 


V 1416 + 4446 


0.400 


541 


31.0 


I-V 


73.8 


2130 


MACS J0159.8 - 0849 


0.405 


3265 


17.6 


I-V 


118.0 


8100 


MACS J2228.5 + 2036 


0.412 


3285 


20 


I-V 


137.7 


6070 


MS 0302.7 + 1658 


0.424 


525 


10.0 


I-V 


59.0 


635 


MS 1621.5 + 2640 


0.426 


546 


30.0 


I-F 


118.0 


3280 


MACS J0417.5 - 1154 


0.440 


3270 


12 


I-V 


138.0 


7400 


MACS J1206.2 - 0847 


0.440 


3277 


23 


I-V 


138.0 


11720 


RX J1347.5 - 1145 


0.451 


3592 


57.5 


I-V 


128.0 


62700 


V 1701 +6414 


0.453 


547 


49.0 


I-V 


64.0 


2745 


CL 1641 +4001 


0.464 


3575 


45.0 


I-V 


49.0 


1040 


MACS J1621.4 + 3810 


0.465 


3254 


9.7 


I-V 


78.0 


1600 


MACS J 1824. 3 +4309 


0.487 


3255 


14.8 


I-V 


84.0 


530 


MACS J1311.0- 0311 


0.492 


3258 


14.8 


I-V 


79.0 


2100 


V 1525 + 0958 


0.516 


1664 


50 


I-V 


79.0 


2100 


MS 0451.6 - 0305 


0.539 


529,902 


56 


I/S-V 


98.4 


16850 


MS 0015.9 + 1609 


0.541 


520 


67.0 


I-V 


98.4 


16200 


MACS Jl 149.5 + 2223 


0.544 


1656,3589 


38 


I-V 


148.0 


9400 


MACS J 1423. 8 + 2404 


0.545 


1657 


18.5 


I-V 


79.0 


3600 


MACS J0717.5 + 3745 


0.548 


1655,4200 


78 


I-V 


144.0 


29000 


V 1121 +2327 


0.562 


1660 


70.0 


I-V 


69.0 


2050 


SC 1120- 1202 


0.562 


3235 


68 


I-V 


49.0 


730 


RX J0848.7 + 4456 


0.570 


927,1708 


184.5 


I-V 


30.0 


850 


MACS J2129.4- 0741 


0.570 


3199 


17.6 


I-V 


98.0 


3000 


MS 2053.7 - 0449 


0.583 


551,1667 


88 


I-V 


54.1 


2150 


MACS J0647.7 + 7015 


0.584 


3196 


19.2 


I-V 


88.5 


3170 


RX J0956.0 + 4107 


0.587 


5294 


17.2 


I-V 


64.0 


500 


CL 0542.8 -4100 


0.634 


914 


50 


I-F 


78.7 


2220 


RCS J1419.2 + 5326 


0.640 


3240 


9.7 


S-V 


44.0 


470 


MACS J0744.9 + 3927 


0.686 


3197,3585 


40 


I-V 


98.0 


6100 


RX J1221.4 + 4918 


0.700 


1662 


78 


I-V 


78.7 


2900 


RX J1113.1 - 2615 


0.730 


915 


103 


I-F 


39.4 


1200 


RX J2302.8 + 0844 


0.734 


918 


108 


I-F 


54.0 


1600 


MS 1137.5 + 6624 


0.782 


536 


117 


I-V 


49.2 


4150 


RX J1317.4 + 2911 


0.805 


2228 


110.5 


I-V 


24.5 


240 


RX J1350.0 + 6007 


0.810 


2229 


58 


I-V 


64.0 


750 


RX J1716.4 + 6708 


0.813 


548 


51 


I-F 


54.0 


1520 


RX J0152.7 - 1357 S 


0.828 


913 


36 


I-F 


52.7 


570 


MS 1054.4 - 0321 


0.832 


512 


80 


S-F 


78.7 


10000 


RX J0152.7 - 1357 N 


0.835 


913 


36 


I-F 


58.0 


830 


1WGAJ1226.9 + 3332 


0.890 


932,3180 


9.5 


S-V 


64.0 


2400 


CL 1415.1 +3612 


1.030 


4163 


89 


I-V 


39.4 


1320 


RDCS J0910 + 5422 


1.106 


2227,2452 


170 


I-V 


24.6 


440 


RXJ1053.7 + 5735 E* 


1.134 


4936 


94 


S-V 


28.2 


300 


RX J1053.7 + 5735 W* 


1.134 


4936 


94 


S-V 


28.2 


450 


RDCS J 1252 -2927' 


1.235 


4198,4403 


188.4 


I-V 


34.5 


850 


RDCS J0849 + 4452* 


1.261 


927,1708 


184.5 


I-V 


23.6 


360 


RDCS J0848 + 4453 


1.273 


927,1708 


184.5 


I-V 


19.7 


130 



Notes: a observation identification number; b effective exposure time after removal of high background intervals; c detector (ACIS-I or -S) and 
telemetry (FAINT or VFAINT) used; d extraction radius; c number of net detected counts in the 0.3 - 10 keV band; * denotes clusters for 
which we also use XMM-Newton observations (see Table|2j. 
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CTI correctiorQ to the observations taken when the tempera- 
ture of the focal plane was 153 K. This procedure allows us to 
recover the original spectral resolution that is partially lost be- 
cause of the CTI. The correction applies only to ACIS-I chips, 
since the ACIS-S3 did not suffer from radiation damage. 

For data taken in the FAINT mode we ran the tool 
acis.process.events only to apply the CTI and the time- 
dependent gain correction. From this point on, the reduction 
was similar for both the FAINT and the VFAINT exposures. 
The data were filtered to include only the standard event grades 
0, 2, 3, 4, and 6. We checked visually for hot columns left from 
the standard cleaning. Only in a few cases did hot columns have 
to be removed by hand. We identify the flickering pixels as the 
pixels with more than two events contiguous in time, where 
a single time interval was set to 3.3 s. For exposures taken 
in VFAINT mode, there were practically no flickering pixels 
left after filtering out "bad" events. We finally filtered time 
intervals with high background by performing a 3cr clipping 
of the background level using the script analyze_ltcr\fl 
Removed time intervals always amount to less than 5% of the 
nominal exposure time for ACIS-I chips. Some ACIS-I obser- 
vations show large flares on the ACIS-S3 chip (which is on by 
default but not used in the data analysis), but the corresponding 
time intervals are not removed since the flares do not affect the 
ACIS-I chips. In any case, our spectral analysis is not strongly 
affected by residual flares, since we always compute the back- 
ground from source-free regions around the clusters from the 
same observation (see below), thus taking into account any pos- 
sible spectral distortion of the background itself induced by the 
flares. 

As in Paper I, we performed a spectral analysis extracting 
the spectrum of each source from a region defined in order to 
maximize the signal-to-noise ratio (S/N, see Appendix A. 1 for 
details on how this is computed). This choice of the extrac- 
tion region allows the global properties of the clusters to be 
measured using the majority of the signal. This strategy is opti- 
mized for the highest redshift objects, and it is homogeneously 
adopted for the whole sample. The region of maximum S/N is 
obtained through the following procedure: the spectrum of each 
source is extracted from a circular region around the centroid 
of the photon distribution. For a given radius, we find the center 
of the region that includes the maximum number of net counts 
in the 0.5 - 5 ke V band, where the bulk of the source counts are 
detected. Then, we compute the S/N, repeating this procedure 
for several radii. Finally we choose the extraction radius R ext , 
defined as the radius for which the S/N is maximum. As shown 
in Fig.Q] in most cases R ext is be tween 0.15 and 0.3 times the 
virial radius /?,,„•, estimated, after Evrard et al. (1996), as 



*™" = 3 - 95 ( H^vF ^ )MPC - 



(1) 



Here the virial temperature 7\,,> is approximated with the spec- 
tral temperature T spec measured within R ext , and 



F(z) = (A(z)/A„r I/6 [Q (1 + z? + 1 - O ] 



-1/2 



(2) 
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Fig. 1. Extraction radius (R ex t) versus virial radius (/?,,» ) for the 
whole sample. Lower and upper lines show R ext = 0.15 R V ir and 
R ext = 0.3 R V i r , respectively. 



where A(z) is the density contrast of the virialized halo with 
respect to the critical density. For sim plicity we assume a con - 
stant value here for (A(z)/A ) _1/6 (see lBrvan & Normanll998l) . 
From Fig. Q] we note that the typical extraction region, depend- 
ing both on the redshift through the surface brightness dim- 
ming and on the brightness distribution of each source, does 
not show a clear trend with redshift, with the exception of the 
highest-z bin where R ex , < 0.15 R v i r , The fraction of net counts 
included in the extraction region always amounts to 80 - 90% 
of the total detected for each cluster. We also note that R ext is 
roug hly 3 times the c ore radius measured with a beta model 
(see Ettori et aT]l2003 for the spatial analysis of a subsample of 



4 http://cxc.harvard.edu/ ciao/threads/acisapplycti/ 
3 http:// cxc.harvard.edu/ciao/threads/filter_ltcrv/ 



our clusters). 

For each cluster, we used the events included in each of the 
extraction regions defined above to produce a spectrum (pha) 
file. The background is always obtained from empty regions of 
the chip in which the source is located. This is possible since all 
sources have an extension of less than 3 arcmin, as opposed to 
the 8 arcmin size of the ACIS-I/-S chips. The background file 
is scaled to the source file by the ratio of the geometrical area. 
The background regions should partially overlap with the outer 
virialized regions of the clusters. However, the cluster emission 
from these regions is negligible compared to the instrumental 
background and does not affect our results. Our background 
subtraction procedure, on the other hand, has the advantage 
of providing the best estimate of the background for that spe- 
cific observation. By comparing the count rate in the source 
and in the background at energies higher than 8 keV, we finally 
checked that variations in the background intensity across the 
chip did not affect the background subtraction, where the sig- 
nal from the clusters is null. The response matrices and the 
ancillary response matrices of each spectrum were computed 
respectively with mkacisrmf and mkwarf for the same re- 
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gions from which the spectra were extracted. For those obser- 
vations for which the CTI correction cannot be applied (when 
the temperature of the detector is larger than 153 K), we used 
acisspec instead. 

2.2. XMM-Newton data 

As in Paper I, we used the XMM-Newton data to boost the 
S/N only for the most distant clusters in our current sample, 
namely the clusters at z > 1 . In Table [2] we list the four 
XMM-Newton observations of high redshift (z > 1) clus- 
ters included in our analysis. For each observation we used 
both the European Photon Imaging Camera (EPIC) PN and 
the two MOS detectors. The XMM-Newton observations and 
data reduction relative to RDCS J0849 and the two clumps of 
RX J1053 have already been presented in Paper I, while th ose 
relative to RDCS J1252 are described in lRosati etail |2004). 



3. Results 

In this section we present the main results of our analysis. 
The section is subdivided into five subsections. In Sect. 3.1 
we provide a general description of our spectral analysis and 
a comparison with the previous results obtained in Paper I. In 
Sect. 3.2 we describe the single source analysis and the main 
properties of the sample. In Sect. 3.3 the correlation between 
iron abundance and temperature is discussed. In Sect. 3.4 we 
present our results on the evolution of the average iron abun- 
dance as a function of redshift obtained through two indepen- 
dent methods (combined fits and weighted means). Finally, in 
Sect. 3.5 we present a comparison with the local iron abun- 
dance of the ICM. 



3.1. Spectral analysis 



The spectra were analyzed with XSPEC v 11.3.1 (|Arnaud 
1996 1 ) and fitt e d with a single- temperature mekal model 



Liedahl et al. 



dKaastral 119921 
tween the elements was fixed to 



1995) in which the ratio be- 
the solar value as in 
Anders & Grevessd (119891) . These values for the solar metal- 
licities have more recently been super seded by the new val - 
ues by iGrevesse & Sauval ( 1998 ) and Asplund et al. (2005), 
who introduced a 0.676 and 0.60 times lower iron solar abun- 
dance, respectively (photometric value). However, we prefer 
to report iron abundances in units of solar abundances by 
Anders & Grevessd (11989b since most of the literature still 
refers to t hem. We also perfor med the fits using solar abun- 
dances by I Asplund et al. (2005). The iron abundances in these 
units (reported in the fifth column of Table 3) can be obtained 
with an accuracy of about 10% simply by rescaling t he values 
measured in solar units by lAnders & Grevesse (1989) by a fac- 
tor of 1.6. This shows that we are not affected by the presence 
of metals other tha n iron. Finally, w e model Galactic absorp- 
tion with tbabs JWilms etaljlioooh . 

It has recently been shown that a methylene layer on the 
Chandra mirrors increases the effective area at energies higher 



than 2 keV dMarshall et al.l E004^[ This has a small effect on 
the total measured fluxes, but it may be non-negligible on the 
spectral parameters (i.e., it may artificially reduce the temper- 
atures). In order to correct for it, we introduced a "positive ab- 
sorption edge" (XSPEC model edge) in the fi tting model at 
2.07 keV with t = -0.15 dVikhlinin et ai1 l2005l). 

The fits were performed over the energy range 0.6 - 8.0 
keV. Due to uncertainties in ACIS calibration below 0.6 keV, 
we excluded less energetic photons from the spectral analysis 
in order to avoid systematic bias. The effective cut at high en- 
ergies is generally lower than 7-8 keV, since the S/N for a 
thermal spectrum rapidly decreases above 5 keV. 

The free parameters in our spectral fits are temperature, 
metallicity, and normalization. Local absorption is fixed to the 
Galactic neutral hydrogen co lumn density (Nh in Tab le [3J, 
as obtained from radio data (IDickey & Lockmanl ll990). and 
the redshift to the value measured from optical spectroscopy 
(z in Table |3). We used Cash statistics applied to the source 
plus background], wh ich is preferable for low S/N spectra 
dNousek&Shue1l989l) . 

Eventhough our analysis procedure is similar to the one 
used in Paper I, new calibrations (including the treatment of 
the positive edge at 2.07 keV) may cause some differences in 
the new temperature and iron abundance values compared to 
the previous analysis. To investigate such differences, we show 
in Fig. 12 the temperatures and iron abundances published in 
Paper I, plotted against the new values (presented in this paper) 
for the 17 clus ters observed with Chandra in both samples, plus 
RDCS 1252 dRosati et al.ll2004 . The new best-fit temperatures 
(upper panel) seem to be slightly higher than in Paper I, while 
iron abundances (lower panel) are much less affected by the 
new calibrations. As a further check we recomputed the aver- 
age values of the iron abundance in the same redshift bins used 
in Paper I and we find consistent results (see Fig.0. The only 
noticeable difference is the hint of a decrease with z, well be- 
low the 2cr confidence level in Paper I, which is now slightly 
enhanced. Therefore, the most significant improvements with 
respect to Paper I are due to the new clusters included in the 
current sample. In particular, we recall that the z — 1.2 point 
was entirely dominated by the XMM data on RX J1053, since 
in Paper I there was no statistically significant detection of the 
iron line at z > 1 from Chandra data only. The situation at z > 1 
has improved thanks to the Chandra observations of CL 1415 
and RDCS 1252. Finally, the statistics of the current sample 
mostly improved in the redshift range 0.3 < z < 0.6. 



3.2. Single source analysis 

This section presents the results of the spectral analysis of each 
of the 56 clusters of our sample. The results of the spectral fits, 
referring to the region of radius R exh defined in Sect. 2.1, are 
listed in Table [3] Histograms of the redshift and temperature 



6 http://cxc.harvard.edu/ccw/proceedings/03_proc/presentations/ 
marshall2 

7 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/ 
XSappendixCash.html 
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Table 2. Additional XMM-Newton observations at z > 1 . 



Cluster 


z 


Exp. [ks] a 


Detector b 


Rex, l"T 


Net Counts' 1 


RXJ105346.6 + 573517 E 


1.134 


94.5 


PN+2MOS 


32 


708 


RXJ105346.6 + 573517 W 


1.134 


94.5 


PN+2MOS 


32 


875 


RDCS J1252 - 2927 


1.235 


65.0 


PN+2MOS 


34.5 


1570 


RDCS J0849 + 4452 


1.261 


112.0 


PN+2MOS 


29.5 


630 



Notes: a effective exposure time after removal of high background intervals; 
counts in the 0.3 - 10 keV band. 



detectors used; c extraction radius; number of net detected 




4 6 8 10 

kT (keV) Tozzi et al. (2003) 





0.2 0.4 0.6 

Z FB /Z Q Tozzi et al. (2003) 



Fig. 2. Comparison between the temperature (upper panel) and 
iron abundance (lower panel) values measured in Paper I, and 
those measured in this work after the most recent Chandra cal- 
ibrations have been applied. Dashed lines show the locus of 
equal temperature and abundance values. 



distribution of the sample are shown in Figs. [4] and [5] respec- 
tively. 

The redshift distribution is peaked around z — 0.5, while at 
z > 1 we have only 7 objects (we rec all that we consider the two 
clumps of RX J 105 3 separately, see lHashimoto et al .1 12004 . In 
order to investigate the properties of the sample as a function 



Fig. 3. Average iron abundance in different redshift bins com- 
puted for the same sample of clusters analyzed in Paper I using 
updated calibrations (red circles), compared with the previous 
results (black circles). The plot shows that the new calibrations 
and data reduction have a negligible effect on our results. Only 
clusters with kT > 5 keV are considered. 



of redshift, we divide the clusters in 5 redshift intervals (10 
objects with 0.3 < z < 0.4; 12 objects with 0.4 < z < 0.5; 15 
objects with 0.5 < z < 0.65; 12 objects with 0.65 < z < 1.0, 
and 7 objects with 1.0 < z < 1.3). These intervals (Fig.lH are 
defined in order to have a comparable number of objects in a 
reasonably narrow redshift range. 

As shown in the temperature distribution (Fig. [5j, we sam- 
pled mostly hot clusters (kT > 5 keV), while only 12 are in the 
medium temperature range (3 < kT < 5 keV). We would like 
to point out here that we derived a single spectral temperature 
for the region within R ext . In principle, the spectral temperature 
can be significantly different from the emission-weighted and 
gas-mass-weighted temperature, in the pr esence of a thermally 
complex ICM (e.g. iMazzotta et alj|2004l) . In a few cases, the 
effect of the temperature gradient is strong enough to make the 
single-temperature fit unacceptable. To evaluate the goodness 
of each fit, we computed the x 1 of each best-fit model after 
binning the spectrum to 20 counts per bin. We find that most 
of our spectra are well-fitted by a single-temperature mekal 
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Table 3. Spectral fit results obtained with the tbabs (mekal ) model. 
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Cluster 


z 


kT [keVf 


Z/Z e (And. & Gre.) b 


Z/Z (Aspl.) c N H [cm 


-2 ]d 


X 2 r [d.o.f.] e 


Null3Typ. Prob. f 


MS 1008.1 - 1224 


0.306 




30 +0 0/ 

u,JU -0.06 


47+o.n 

-0.10 


7.26 


X 


10 2 " 


1.23 


[228] 


0.009 


MS 2137.3 - 2353 


0.313 


4.96 ±0.11 


0.35 ± 0.03 


0.56!«« 


3.55 


X 


10 20 


1.23 


;289] 


0.004 


Abell 1995 


0.319 


8.60 ± 0.32 


0.40 ± 0.06 


0.64!» : i» 


1.42 


X 


10 20 


1.22 


343] 


0.003 


MACS J0308.9 + 2645 


0.324 


11.2 + 0.7 


0.37 ± 0.06 


0.70*}° 


1.18 


X 


o 21 


1.07 


;283] 


0.194 


ZwCl 1358.1 +6245 


0.328 


6.70 + 0.26 


0.40 ± 0.06 


0.64+o:jI 


1.92 


X 


10 20 


1.23 


;282] 


0.006 


MACS J0404.6 + 1109 


0.355 


6 Q+° 6 


1 fi +au7 
u.io_ u 


22 +al4 

"•^ -0.20 


1.43 


X 


LO 21 


1.01 


;i65] 


0.451 


RXJ0027.6 + 2616 


0.367 


O 1+2.6 


57 +0 - 27 

—0.19 


1.02!"- 34 


3.86 


X 


10 20 


1.36 


[52] 


0.042 


MACS J1720.2 + 3536 


0.391 


6.30 ± 0.33 


50 

-0.06 


0.83!™ 7 


3.40 


X 


10 20 


0.74 


[189] 


0.997 


ZwCl 0024.0 + 1652 


0.395 


4.38 ±0.27 


7i+ - 20 

"• /J -0.18 


i 99+0.32 
1 - i -0.29 


4.20 


X 


10 20 


1.02 


[128] 


0.410 


V 1416 + 4446 


0.400 


3.50 ±0.18 


0.88 ±0.18 


1 50 +a35 

27 


1.29 


X 


10 20 


1.11 


;98] 


0.219 


MACS J0159.8 - 0849 


0.405 


Q 9+O.6 


0.36 ± 0.05 


0.59^ 


2.08 


X 


q2() 


1.13 


;208] 


0.097 


MACS J2228.5 + 2036 


0.412 


7.9 ±0.6 


"■^-0.07 


0.59!™ 


4.58 


X 


10 20 


0.91 


196] 


0.801 


MS 0302.7 + 1658 


0.424 


A 34+O.56 
H -"-0.44 


40 +a21 


0.58! , 34 


1.11 


X 


o 21 


0.98 


[31] 


0.490 


MS 1621.5 + 2640 


0.426 


6 9 +0 7 


35 +a08 

"• JJ -0.09 


0.63 ±0.14 


3.58 


X 


10 20 


0.90 


[157] 


0.804 


MACS J0417.5 - 1154 


0.440 


11.3 ±0.8 


n 9Q+0.06 
"■^-0.07 


45 +ai ° 


3.86 


X 


10 20 


0.82 


;209] 


0.971 


MACS J1206.2-0847 


0.440 


1 1 o +0 - 7 




25 +ai ° 

" -0.09 


3.72 


X 


10 20 


1.15 


;27i] 


0.045 


RX J1347.5 - 1145 


0.451 


14.0 ± 0.4 


0.33 ± 0.04 


55 +a07 

—0 06 


4.92 


X 


10 20 


1.20 


;436] 


0.003 


V 1701 + 6414 


0.453 


A 77+0-26 
^■""-0.25 


o.54t;i 


°- 85 -a20 


2.46 


X 




0.95 


.123] 


0.649 


CL 1641 + 4001 


0.464 


4.8 ±0.6 


0-48!°;^ 


0.79!^ 


1.10 


X 


10 20 


0.98 


[53] 


0.509 


MACS J1621.4 + 3810 


0.465 


6.5 ± 0.7 


0.13 ±0.08 


0.19+^ 


1.09 


X 


q2() 


1.46 


[67] 


0.549 


MACS J1824.3 +4309 


0.487 


7 2 +22 

'■^-1.3 


0.38!°| 


55 +a31 

"• JJ -0.33 


4.46 


X 


10 20 


1.03 


[33] 


0.423 


MACS J1311.0- 0311 


0.492 


8.0 ±0.9 


0.39 ± 0.09 


°- 60+ 019 


1.87 


X 


10 20 


0.97 


[87] 


0.571 


V 1525 + 0958 


0.516 


c 4+0.4 


0.35 ±0.10 


54 +a21 

U,J -0 20 


2.91 


X 


10 20 


1.06 


107] 


0.315 


MS 0451.6-0305 


0.539 


o 9+0.4 
8.2_ 03 


0.34 ± 0.06 


0.57 ±0.10 


4.97 


X 


10 20 


1.24 


;279] 


0.004 


MS 0015.9 + 1609 


0.541 


q a+0.5 
7 ' -0.3 


r\ "3 q +0.06 
U,J -0.05 


0.53 ± 0.09 


4.07 


X 


qZO 


1.01 


302] 


0.460 


MACS J 1149.5 + 2223 


0.544 


1 2 9 +L2 


2 1 

u,i -0.07 


0.37!»j; 


2.28 


X 


q2() 


1.04 


;254] 


0.317 


MACS J1423.8 + 2404 


0.545 


7.3_ a5 


0-31^08 




2.38 


X 


2 " 


1.11 


;i35] 


0.174 


MACS J0717.5 + 3745 


0.548 


13.3 ±0.7 


0-18!^ 


29 +0 " 

"""-0.07 


7.04 


X 


10 20 


1.14 


;404] 


0.029 


V 1121 + 2327 


0.562 


5.2 ± 0.5 


.-. — i A 1 A 

27 +,,1(l 

"■^■'-0.08 


43+o.i5 

-0 07 


1.32 


X 


10 20 


0.99 


;i02] 


0.499 


SC 1120- 1202 


0.562 


5 7 + '- 3 


^* i A OA 

23 

"•"-0.17 


33 +o:33 

—0 21 


5.19 


X 


10 20 


1.13 


[45] 


0.256 


RX J0848.7 + 4456 


0.570 


3.4 ±0.30 


51 +a23 

"• J -0.19 


79+ 36 

u - -0.30 


2.63 


X 


10 20 


1.33 


[46] 


0.068 


MACS J2129.4 - 0741 


0.570 


o -7+0.7 


o sr - 08 

" J '-0.11 


0.82 ±0.14 


4.82 


X 


10 20 


0.85 


332] 


0.890 


MS 2053.7 - 0449 


0.583 


5 - 68 -a47 


16 +0 " 

"• lu -0.10 


24 +ai8 

-0.16 


5.02 


X 


10 20 


1.03 


300] 


0.402 


MACS J0647.7 + 7015 


0.584 


15.51 2 - 7 


< 0.10 


<0.15 


5.64 


X 


o 2 " 


0.91 


.118] 


0.737 


RX J0956.0 + 4107 


0.587 


7 4 2 5 

' '^-1.45 


1 3+ u -28 
"■ 1J -0.13 


21 +a45 

"•^ -0.21 


1.14 


X 


q2() 


1.08 


[27] 


0.353 


CL 0542.8 - 4100 


0.634 


7 9 +u 

'• 7 -0.8 


20 +012 

u,i -0.09 


0-3 it!? 


3.73 


X 


10 20 


1.19 


313] 


0.081 


RCS J1419.2 + 5326 


0.640 


4 1 +0.7 


1 i+O- 21 


0.16^i 


1.18 


X 


o 20 


1.35 


[21] 


0.131 


MACS J0744.9 + 3927 


0.686 


9.2 ±0.6 


0.28 ± 0.06 


0.46!» : J» 


5.71 


X 


10 20 


1.22 


[193] 


0.021 


RX J1221.4 + 4918 


0.700 


8 4 +0 9 


29 +0 B 


°- 48+ 019 


1.47 


X 


10 20 


0.99 


341] 


0.534 


RX Jl 1 f 3.1 -2615 


0.730 




. — [ A OA 

43+0.20 

"■^ -0.17 


0.68t» 


5.50 


X 


10 20 


0.69 


[59] 


0.966 


RX J2302.8 + 0844 


0.734 


s.otiJ 


12 +a ' 6 

"■ lz --0.12 


19 - 26 

" -0.19 


4.85 


X 


10 20 


0.85 


[81] 


0.833 


MS 1137.5 + 6624 


0.782 


6.8 ± 0.5 


0.26 ±0.11 


43+o.is 

"•^ -0.17 


1.21 


X 


20 


1.09 


[151] 


0.204 


RX J1317.4 + 2911 


0.805 


4 5 +L4 


3'5+°- 46 

"• JJ -0.30 


49+0.67 
"•^ -0.46 


1.10 


X 


10 20 


0.71 


[14] 


0.768 


RX J1350.0 + 6007 


0.810 


4 4+0-7 
H - H -0.6 


56 +0 27 

U - J -0.18 


96 +a35 

v - y -0.33 


1.80 


X 


10 20 


1.26 


[51] 


0.103 


RX J1716.4 + 6708 


0.813 


6.9!™ 


49+°- 18 

"•^ 7 -0.16 


79+ - 29 

' 7 -0.26 


3.72 


X 


10 20 


0.77 


[76] 


0.934 


RX J0152.7 - 1357 S 


0.828 


8- 7 ^8 


< 0.22 


<0.37 


1.54 


X 


10 20 


1.06 


[32] 


0.376 


MS 1054.4-0321 


0.832 


7 5 +a7 

'■-'-0.4 


23 +a07 

" -0.08 


0.38!°! 2 


3.61 


X 


10 20 


1.15 


338] 


0.056 


RX J0152.7 - 1357 N 


0.835 


6 7 +12 

"•'-1.0 


17+ 019 

-0.16 


0.25!°i? 


1.54 


X 


10 20 


0.87 


[45] 


0.715 


1WGAJ1226.9 + 3332 


0.890 


12.9 L4 2 


02 +012 


< 0.21 


1.38 


X 


20 


1.06 


[97] 


0.335 


CL 1415.1+3612 


1.030 


7 +0 - 8 ' 


24 +015 

"■^-0.13 


40 +0 - 23 

u - w -0.21 


1.09 


X 


10 20 


0.79 


;63] 


0.891 


RDCS J0910 + 5422 


1.106 


6 4 +1 5 

-1.2 


<0.14 


< 0.21 


2.10 


X 


10 20 


0.94 


[28] 


0.559 


RXJ1053.7 + 5735E 


1.134 


3 4 +05 


51 +014 

-0.15 


°- 58 -a22 


5.7 x 10 19 


1.11 


[29] 


0.317 


RXJ1053.7 + 5735 W 


1.134 


7 2 +l2 

' -0.9 


32 +014 

"-' z --0.15 


0.51!°; 23 


5.7 x 10 19 


0.55 


[31] 


0.981 


RDCS J 1252 - 2927 


1.235 


7 9+0.4 
'• -0.6 


35+ 06 

U - JJ -0.09 




5.95 


X 


20 


1.01 


[55] 


0.454 


RDCS J0849 + 4452 


1.261 


6 2 +1 


16+ 013 

u - lu -0.14 


021 +0.27 
"•^• 1 -0.19 


2.63 


X 


10 20 


0.51 


[24] 


0.977 


RDCS J0848 + 4453 


1.273 


3 4 +25 


22 +L33 

"•^ -0.22 


22 +L38 

"•^ -0.22 


2.63 


X 


10 20 


1.04 


[11] 


0.411 



Notes: a temperature; iron abundance in solar units bv lAnders&Grevessd dl989h and c by I Asplund etai] J2005b : d local column density, 
always fixed to the Galactic value by Dickey & Lockman (1990); c reduced chi-square and degrees of freedom obtained after binning the 
spectra to 20 counts per bin; f null-hypothesis probability. Errors refer to the lcr confidence level. 
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Fig. 4. Histogram of the redshift distribution of the sample. The 
solid line refers to the whole sample, while the dotted line dis- 
plays only clusters with kT < 5 keV (see TableO. The vertical 
dotted lines indicate the five redshift intervals selected for the 
combined spectral analysis, as described in the text. 



Fig. 6. Null-hypothesis probability for the single-temperature 
best fits as a function of net detected counts for the whole sam- 
ple. 



10 



II 



10 

kT (keV) 



Fig. 5. Histogram of the temperature distribution of the sam- 
ple. The vertical dotted line separates the moderate temperature 
(kT < 5 keV) clusters from the high temperature [kT > 5 keV) 
clusters. 



model (see 7th column in Table 0. However, when the num- 
ber of the net detected counts becomes larger than ~ 10 4 , the 
quality of the fits drops dramatically for about half of the clus- 
ters (see Fig. [6]). If we consider a 1% null-hypothesis probabil- 
ity as the threshold for an acceptable fit, then we must reject 
the single-temperature model within R ext for 5 clusters in our 
sample (namely MS 2137, A 1995, ZW 1358, RX J1347, and 



MS 045 1 ). Note that we can also fit clusters with very disturbed 
morphology and high S/N (e.g. MACS J0717) with a single- 
temperature model, due to the very high temperatures involved, 
which provide composite spectra with much less features than 
spectra with low-temperature components. Indeed, the single- 
temperature model fails mostly when a strong cool-core with 



temp eratures lower than 3 keV is present (e.g. Maz zotta et al 
2004), while it is still acceptable if the temperature range is 
well above this threshold. 

Since we are focusing here on metallicity, we made a closer 
investigation of the best-fit Zp e values for the clusters with the 
lowest null-hypothesis probability. Since these clusters are also 
the ones with the highest S/N, we were able to perform a spa- 
tially resolved spectral analysis for about four concentric an- 
nuli. We find that the best-fit Z Fe value measured with a sin- 
gle temperature mekal model within R ext , is representative of 
the inner 400 kpc and is not dominated by the central bin. 
This result, implying that the presence of temperature gradient 
does not dramatically affect the measurements of iron abun- 
dance, is reinforced by the spectral simulations described in the 
Appendix. Therefore we used the single-temperature best-fit 
values for all the clusters in our sample. The attempt to model 
the evolution of Zp e separately in the inner 100 kpc and the 
outer regions is deferred to a future paper. 

Finally, we show in Fig.|7]the distribution of temperatures 
in our sample as a function of redshifts (error bars are at the 
lcr confidence level). The Spearman test shows no correlation 
between temperature and redshift (Spearman's rank coefficient 
of r s = -0.095 for 54 degrees of freedom, probability of null 
correlation p = 0.48). Fig. [7] shows that the range of tempera- 
tures in each redshift bin is about 6-7 keV. Therefore, we are 
sampling a population of medium-hot clusters uniformly with 
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Fig. 7. Temperature plotted vs redshift for the whole sample. 
Shaded areas show the rms dispersion around the weighted 
mean. 

redshift, with the hottest clusters preferentially in the redshift 
bin0.4<z<0.6. 

The relations between temperature and iron abundance at 
different redshifts are shown in Fig. [8] For three clusters we can 
only derive upper limits on the iron abundance of the ICM, two 
of them at z < 0.8. For five clusters we measure positive iron 
abundances, which are still consistent with no detection at the 
lcr c.l.. Overall, we detect the presence of the iron line in the 
large majority of the clusters, and measure the iron abundance 
with a typical error of 30% at the lcr c.l. for z < 0.6, and 50% 
or larger at z > 0.6. 

3.3. The iron abundance-temperature correlation 

Our analysis (Fig.|9]l suggests higher iron abundances at lower 
temperatures in all the redshift bins. This trend is somewhat 
blurred by the large scatter. We find a more than 2<x negative 
correlation for the whole sample, with a Spearman's rank co- 
efficient of r s = -0.31 for 54 degrees of freedom (probability 
of no correlation p = 0.018). The correlation is more evident 
when we compute the weighted average of the metallicity in six 
temperature intervals (see Table|4|, as shown by the shaded ar- 
eas in Fig. [9] The weighted mean is computed as T wa = ST/w, 
in each temperature bin, where w; = 1 jcrfj '1,(1 /erf), and cr ; is 
the lcr error on the single measurement. We find that in each 
bin the scatter of the best-fit values around the mean is com- 
parable to the statistical errors on the single measurements (re- 
duced xl — 1 assuming a constant Zp e in the bin), with the 
exception of the third and sixth bins, where the intrinsic scatter 
is larger (xl - 3)- 

This trend is similar to what is found in the A SCA 
data of nearby c l usters by Baum gartner et al. 20051 (see 
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□ 0.65<z<1.0 - 
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also Arnaudet al] 1992 ; Mushotzkv & Loewensteinl 1 1997 



Finoguenov et alj 2001 ). In their paper, the average iron abun 



Fig. 8. Iron abundance-temperature plots for the whole sam- 
ple. The four panels show each redshift bin separately (last two 
bins in the fourth panel). The dashed line represents the best-fit 
metallicity-temperature relation (Z/Z 0.88 T~ 0A1 ) referring 
to the whole sample. Error bars refer to lcr confidence level. 
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Fig. 9. Scatter plot of best-fit iron abundance values (with- 
out error bars) versus temperature for the whole sample. The 
dashed line represents the best-fit metallicity-temperature re- 
lation (Z/Z Q =i 0.88 T~ 0A1 ). Shaded areas show the weighted 
mean (blue) and average iron abundance with rms dispersion 
(cyan) in 6 temperature bins (see TablelU. 



Fig. 10. Best-fit confidence contour plot for the Z — T relation 
modelled with a power law of the form Z(T) = Zf e (0) T~ C * T . 
Inner contours display the 1, 2, and 3<x c.l. using errors on the 
weighted mean, while outer thick contours are obtained using 
the rms dispersion. 



Table 4. Average iron abundance calculated in different tem- 
perature bins. 



kT [keV] a 


z/z Q b 


AZ/Z Q C 




(weighted mean) 


(rms) 


3.9 [10] 


0.542 ± 0.060 


0.23 


5.4 [8] 


0.329 ± 0.024 


0.10 


6.7 [11] 


0.322 ± 0.027 


0.17 


7.5 [10] 


0.317 ±0.030 


0.11 


8.8 [9] 


0.350 ± 0.024 


0.16 


12.7 [8] 


0.244 ±0.021 


0.14 



Notes: a average temperature in each bin (the number of clusters in 
each bin is shown in parenthesis); b weighted mean of the iron 
abundance; c rms dispersion. 



dance for kT > 5 keV is constant and equal to Z ^ 0.3 Z Q , while 
it rises to very high values (well above 0.4 Z Q ) in the tempera- 
ture range 2-3 keV, and drops below 0.3 Z for kT < 2 keV. 
It is worth noting, however, that this behavior is different from 
that of Ni and a-elements. Here, we confirm that in our high- 
redshift sample the average measured iron abundance rises be- 
low 5 keV. A simple power law fit of the form 



Z(T) =Z Fe (0)T- aT 
to our data gives Zf e (0) 

2 



0.88^ Z and a T 



(3) 



-0.47+01 with 



Xy — 1.8 for 5 degrees of freedom (see the confidence con- 



tours in the slope-normalization plane in Fig. [TO). This trend is 
consistent with what is found from ASCA data of a local sam- 



ple of cluster s for kT > 3 keV (see Baumgartne r et al.ll2005 



perature and metallicity is observed when a spatially-resolved 
analysis of the ICM can be performed in a singl e object, as in 



the case of the core of the Perseus cluster (see Sand ers et al 



Hornerll2005l) . Interestingly, a similar correlation between tem- 



2004). However a physical explanation for this behaviour is 
still missing. 

We recall that the measure of the iron abundance in local 
clusters is based on both the K-shell and the L-shell complex, 
at 6.7 - 6.9 keV and between 1 and 2 keV, respectively. It has 
been pointed out that a diagnostic based mostly on the L-shell, 
as in the case of spectra with a significant l ow temperature com- 
ponent, is more uncertain (IRenzinill 19971) . In our high redshift 
sample, we expect to be sensitive only to the iron K-shell com- 
plex. Indeed, when we fit the spectra cutting energies below 
2 keV rest-frame, we find that the best-fit metallicities are con- 
sistent with those found within the 0.6 - 8 keV range (observed 
frame) used throughout the paper, as shown in Fig.Ql] We no- 
tice that when the low energy range is removed, errors on the 
temperatures become larger, with a clear tendency to higher 
values. This trend is expected since higher temperature spec- 
tral shapes can be accommodated more easily than lower ones. 
On the other hand, iron abundances hardly change, except for 
the two most iron-rich clusters, ZW 0024 and V 1416. This may 
indicate that in these two objects the high iron abundance could 
be associated with a low temperature component. However, in 
these cases a separate analysis in two annuli is possible and we 
do not find a clear enhancement of the iron abundance in the 
inner regions. 

The observed trend between Zp e and T is still a matter of 
debate. It may be linked to the observed decrease in the star 
formation efficiency with increasing cluster mass as reported 
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Fig. 11. Best-fit temperatures (upper panel) and iron abun- 
dances (lower panel) obtained using the 1.5-8 keV energy 
range compared to the values obtained from the 0.6 - 8 keV 
(this paper). Dashed lines show the locus of equal temperatures 
and abundances. Only clusters with more than 2000 net counts 
(see Table[T|i are considered here. 



in lLin et a"L"l(l2003l) . The modelization of this trend goes beyond 
the scope of this paper. 

3.4. The evolution of the iron abundance via 
combined spectral analysis 

In Fig.[T2]we show the iron abundance best-fit values for all the 
sources in the sample. When focusing on the highest redshift 
bin, it is worth noticing that at z > 1 we find clear detections of 
the Fe K line in the spectra of CL J1415 (z = 1.030, Z Fe > 
at the 90% c.L), of RDCS J1252 (z = 1.235, Z Fe > 0.2 Z at 
the 90% c.L), and of the two clumps of RX J1053 (z = 1.134, 
Zp e > 0.1 Z Q at the 90% c.L). In the spectra of the four other 
clusters, we do not have separate detections of the iron line, but 
all measurements are consistent within lcr with Z Fe 0.3 Z . 
Therefore, at present, we have a much better estimate of the 
metal content of clusters at z — 1 than in Paper I, where the 
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Fig. 12. Iron abundance plotted versus redshift for the 56 clus- 
ters of the sample. Shaded areas show the rms dispersion 
around the weighted mean of the iron abundance in the 5 red- 
shift bins defined in the text. Error bars refer to a 1 cr confidence 
level. 

iron line at z > 1 was only firmly detected in the two clumps of 
RXJ1053. 

We find a ~ 3<x negative correlation between iron abun- 
dance and redshift, with a Spearman's rank coefficient of r s = 
-0.40 for 54 degrees of freedom (probability of a null corre- 
lation p = 0.0023). This correlation is stronger than the weak 
hint (less than a 2<x c.L) of anticorrelation found in Paper I. We 
verified that, if we calculate the Spearman's rank coefficient 
only for the 19 clusters analyzed in Paper I, using the newly 
revised temperature and abundances (consistent with the old 
ones as shown in Fig. |2j, we again obtain very weak correla- 
tion (r s = -0.25 for 17 degrees of freedom, probability of no 
correlation p = 0.30), which is therefore consistent with the 
results reported in Paper I. This confirms that our new results 
should not be ascribed to updated calibrations, but rather to the 
larger size of the sample, particularly at z < 0.5. 

The decrease in Z Fe with redshift becomes more evident 
by computing the average iron abundance as determined by a 
combined spectral fit in a given redshift bin. This technique is 
similar to the stacking analysis often performed in optical spec- 
troscopy, where spectra from a homogeneous class of sources 
are averaged together to boost the S/N, thus allowing the study 
of otherwise undetected features. In our case, different X-ray 
spectra cannot be stacked due to their different shape (differ- 
ent temperatures). Therefore, a simultaneous spectral fit is per- 
formed leaving temperature and normalization free to vary for 
each object, and using a unique metallicity value for all the 
clusters in the subsample. 

We note that the scatter of the best-fit values around the 
mean in each redshift bin is, in some cases, larger than the typ- 
ical statistical errors on single measurements. This is expected 
on the basis of the Z - T correlation, as found in Sect. 3.3. The 
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Table 5. Average iron abundance in different redshift bins re- 
sulting from combined fit and weighted mean. 





z/z b 


z/z c 




AZ/Z Q d 




(combined fit) 


(weighted mean) 


(rms) 


0.206 [9] 


n 477+0.003 

"•^'-0.011 


0.416 


± 0.009 


0.08 


0.350 [10] 


0.387^™ 


0.379 


±0.019 


0.20 


0.447 [12] 


0.33<£$™ 


0.318 


± 0.020 


0.11 


0.572 [15] 


306+ 017 

u.juu 0033 


0.260 


+ 0.020 


0.13 


0.787 [12] 


0.244 ± 0.025 


0.251 


± 0.035 


0.15 


1.167 [7] 


265+ 05 


0.28 + 


0.048 


0.15 



Notes: a average redshift of each bin (the number of clusters in each 
bin is shown in parenthesis); b iron abundance from combined 



fit with lcr errors; 1 
dispersion. 



iron abundance from weighted mean; rms 



reduced x 2 obtained by assuming that measurements are scat- 
tered around the weighted average is between 2 and 3 in the first 
two bins, implying the presence of intrinsic scatter comparable 
to the typical statistical error, while it is ~ 1 above z — 0.5, 
since here the typical statistical error is larger than the intrinsic 
scatter component. By assuming a unique value of Z Fe in the 
combined fit, however, we intend to provide an average value 
of the metallicity over large cluster volumes as a function of 
redshift. 

A plot of the combined iron abundance measured in each 
redshift bin is shown in Fig. Qj] To verify the robustness of 
our results, we computed the weighted average from the sin- 
gle source fits in each redshift bin. The best-fit values re- 
sulting from the combined fits are always consistent with the 
weighted means (listed in Table |5]l within lcr, except for the 
bin at z ~ 0.6, which is lower. We also checked that, if the two 
clusters with the highest Z Fe (i.e. ZW 0024 and V1416, see 
Sect. 3.3) are removed from the fit, the average Z Fe value in the 
first redshift bin is only slightly changed (at the level of ~ 3%). 

From a visual inspection of Fig. [13] we notice a continuous 
trend of decreasing iron abundance from z ^ 0.3 to z ^ 1.2. 
While a constant value (Zf e ) ^ 0.25 Z is a good fit at z > 
0.5, the iron abundance is significantly higher at z < 0.5, the 
redshift range over which the statistics of our sample increased 
most with respect to Paper I. In addition, we now have a firm 
measurement of the average iron abundance at redshift z — 0.8 
and z — 1 .2, reinforcing the results of Paper I; in particular, the 
iron abundance in the most distant clusters is still consistent 
with the value Zf e = 0.3 Z Q within lcr. 

Given the negative correlation between iron abundance and 
temperature found in our sample (see Sect. 3.3), we first veri- 
fied whether the evolution with redshift is due to the presence 
of low-temperature clusters in the low redshift bins. By exclud- 
ing clusters with kT < 5 keV from the plot, the same results are 
obtained, as expected from the lack of correlation betwen the 
average temperature and redshift found in Sect. 3.2. 

We also note that this trend points towards (Z Fe ) ~ 
0.5 Z Q at low-z, which is higher than the often-reported value 
Zp e ^ 0.3 Z Q . The reason is that we are measuring Z Fe in 
the inner regions of clusters, where it reaches values signifi- 
cantly higher than 0.3 Z , particularly in cool-core clusters (see 




Fig. 13. Mean iron abundance from combined fits within five 
redshift bins defined in the text (red circles) compared with the 
weighted average of single-source measurements in the same 
bins (black squares). The triangles at z — 0.2 are based on the 
low-z sample described in Sect. 3.5. Error bars refer to the lcr 
confidence level. Shaded areas show the rms dispersion. The 
dashed line indicates the best fit over the 6 redshift bins for a 
simple power law of the form (Zf e ) = Z Fe (Q)(l + zT" z with 
a 7 ^ 1.25. 



Vik hlinin et al.ll2005l) . which constitute about 2/3 of the local 
population. In order to show that the average iron abundance 
in low-z clusters, when analyzed with our procedure, confirm 
the trend seen in the high-z sample, we add a point at (z) - 0.2 
including 9 clusters, as described in detail in Sect. 3.5. 

A fit with a constant iron abundance value over the entire 
redshift range is unacceptable. If we model the evolution (in- 
cluding the additional low-z point, using the values from the 
combined fits) with a simple power law as 



(Z Fe )^Z Fe (Q)(l + z y 



(4) 



the best-fit values obtained are Z Fe (Q) = 0.54 + 0.04 Z and 
a z = 1.25 ± 0.15 (xl = 0.9), implying a decrease by a factor of 
~ 2 between z — 0.3 and 1.2. The evolution (a z < 0) is signif- 
icant only at lcr level when the rms dispersion is used instead 
of the errors on the combined fits (see confidence contours in 
Fig.fPfli. However, the rms dispersion is greatly overestimating 
the uncertainties on the average values. Consistent best-fit re- 
sults are obtained if the lowest redshift point is excluded from 
the fit or if weighted mean values are used. 

3.5. The "local" iron abundance of the ICM 

An important issue to address is how our findings at z <: 0.3 
compare with the local iron abu ndance, which is tradition- 
ally quoted to be Z Fe 0.3 Z Q dRenzini 1997 ) in the units 
following Anders & Grevesse (1989). The spatially-resolved 
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Fig. 14. Best-fit confidence contours plot for the metallicity as a 
function of redshift (including the additional low-z point) mod- 
elled with a power law of the form Z(z) = Zf e (0)(l + z) _Q; *. 
Thick inner contours display the 1, 2, and 3cr c.l. using the 
results of the combined fits, while thin outer contours are ob- 
tained using the rms dispersion. 



spectroscopy of local clusters obtained with the Chandra and 
XMM-Newton satellites shows a complex distribution of the 
metals within the inner regions. The large differences in iron 
abundances and gradients from cluster to cluster lessens the 
meaning of the adoption of a single canonical value for the 
average iron content of the ICM in the local Universe. The 
value Zf e 0.3 Z is typically observed only in non cool- 
core clusters or i n the outer regions (> 100 kpc) of cool- 
core clusters (s ee Tamura et al.l 120041 for XMM-Newton and 
Vikhlinin et al. 20051 for Chandra data). The central peak of 
iron abundance reaches values of Zf e ^ (0.6 - 0.8) Z in the 



core s of cool-core clusters, which have a typical s ize of 100 kpc 
(see lDe Grandi et al.ll2004t I Vikhlinin et alj|2005h . whereas Z Fe 



decreases to ^ 0.3 Z in the outer regions. On the other hand, 
the iron abundance appears to be constant, Z Fe ^ 0.2 - 0.3 Z , 
in non cool-core clusters. As a result, particular care should be 
used when comparing our measurements with the local values 
of Zf e from the literature. 

Since the extrapolation of the average Zp e at low-z points 
towards Zf e (0) ^ 0.5 Z , we need to explain the apparent dis- 
crepancy with the oft-quoted canonical value (Zf e ) ^ 0.3 Z . 
As mentioned in Sect. 3.4, the discrepancy is due to the fact that 
our average values are computed within r ^ 0.15 7?, ,, , where 
the iron abundance is boosted by the presence of metallicity 
peaks often associated to cool cores. The regions chosen for 
our spectral analysis, are larger than the typical size of the cool 
cores, but smaller than the typical regions adopted in studies of 
local samples. 

This can be proved by analyzing the inner regions (r < 
0. 15 R V i r ) of a sample of clusters at z < 0.3. We selected a small 



subsample of 9 clusters at redshift 0.1 < z < 0.3, including 7 
cool-core and 2 non cool-core clusters, a mix that is represen- 
tative of the low-z population. These clusters, listed in Tabled 
are presently being analyzed for a separate project aimed at ob- 
taining spatially -resolved spectroscopy (Baldi et al., in prepa- 
ration). Here we analyze a region within r = 0.15 R v t r in order 
to probe the same regions probed at high redshift. We used this 
small control sample to add a low -redshift point in our Fig.[T3l 
which extends the Zf e evolutionary trend. 



4. Discussion 

The main result of this work is that the cosmic average of Zf e 
in the central regions (R < 0.3 /? V!r ) of clusters significantly de- 
creases with redshift out to z — 0.5, remaining constant out to 
z — 1.3 at the level of Zf e ^ 0.25. Given the complex ther- 
mal and chemical structures observed in bright local clusters, a 
main concern is whether our analysis might be affected by evo- 
lution in the occurrence of temperature/metallicity gradients in 
the cluster population. The assumption of a single-temperature 
mekal model for the inner 0.3i? Vi > may well be too simplis- 
tic, and it may introduce systematic biases in the recovered Zf e 
values. In order to clarify this issue, we investigate possible bi- 
ases in our fitting procedure in the Appendix using a large set 
of simulated spectra in the typical S/N regime of our high-z 
clusters. We find that different S/N values do not introduce any 
significant bias. In particular, spectra with lower S/N (occurring 
mostly at high-redshift) tend to give slightly higher Zf e com- 
pared to the input values, therefore opposite to the observed 
trend (see Fig. A.l). 

In the Appendix we also investigate the cases of a two- 
temperature ICM with a single Zp e and of a two-temperature 
ICM with higher Zf e associated with the colder component, 
analyzed with a single-temperature mekal model. Here the 
key parameter is the ratio of the emission measure of the two 
components. This quantity is difficult to model, so that here we 
assume a few representative values ranging from 0.3 to 0.75. 

In both cases we measured slightly higher temperatures at 
higher redshifts, due to the fact that for high-redshift clusters, 
the signature of the cold component is partially redshifted be- 
low the adopted energy range (E > 0.6 keV). We also find a 
mild trend toward lower Zf e at z ~ 1 compared to z ~ 0.6. 
This effect is limited to be < 30%, and therefore it cannot fully 
explain the observed decrease even under the extreme assump- 
tion that all the clusters in the sample had steep temperature 
gradients with a central abundance peak. To sum up, these sim- 
ulations show that the presence of temperature and abundance 
gradients, // their occurrence is constant through the popula- 
tion of clusters at different redshifts, does not introduce signif- 
icant bias into our measure of the evolution of Zf e . 

Furthermore, we investigate whether the evolution of Zf e 
could be due to an evolving fraction of clusters with cool cores, 
which are known to be associated with iron-rich cores (see 
De Grandi et al. 2004) and which amount to more than 2/3 of 



the local clusters (see Ba uer et al ]l2005h . For example, the evo- 
lution of the mass in iron in the central peak, which is about 
20 - 30% of the total, may be associated with the star forma- 
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Table 6. Spectral fit results for the low-z sample with the tbabs (mekal ) model. 



Cluster 


z 


kT [keVf 


z/z Q b 


N H [cm- 2 Y 


XI- [d-o.f.] d 


Null Hyp. Prob. e 


Abell 1413 


0.143 


7.1 ±0.1 


0.39 ± 0.02 


2.18 x 10 20 


1.55 [477] 


io- 13 


Abell 907 


0.153 


5.2 ±0.1 


0.53 ± 0.03 


5.36 x 10 20 


1.37 [431] 


io- 7 


Abell 2104 


0.155 


13.9 ±0.5 


53 +0 04 


8.69 x 10 20 


1.35 [423] 


10~ 6 


Abell 2218 


0.176 


7.9 ±0.3 


0.26 ± 0.03 


3.26 x 10 20 


1.09 [347] 


0.111 


Abell 963 


0.206 


7.0 ± 0.2 


0.43 ± 0.03 


1.40 x 10 20 


1.09 [341] 


0.113 


Abell 2261 


0.224 


7 c+0.4 
' '-'-0.2 


5 j +0.03 

"""-0.05 


3.28 x 10 20 


1.07 [329] 


0.169 


Abell 2390 


0.228 


9.1 ±0.1 


40 +0 03 

u '^ u -0.02 


6.81 x 10 20 


1.50 [479] 


io- 12 


Abell 1835 


0.253 


7.2 ±0.2 


Q 41 +0.03 


2.32 x 10 20 


0.99 [291] 


0.527 


ZwCl 1021.0 + 0426 


0.291 


6.2 ±0.1 


39+ 03 

"■ J -0.03 


3.02 x 10 20 


1.51 [396] 


io- 10 



Notes: a temperature; iron abundance in solar units bv lAnders & Grevessd dl989h : c local column density, always fixed to the Galactic value 
by Dickey & Lockman (1990); d reduced chi-square and degrees of freedom obtained after binning the spectra to 20 counts per bin; c 
null-hypothesis probability. Errors refer to the lcr confidence level. 



tion p roduct of the central galaxy alone (see iDe Grandi et al 



2004 



In order to use a simple characterization of cool-core clus- 
ters in our high-z sample, we computed the ratio of the fluxes 
emitted within 50 and 500 kpc (C = fir < 50 kpc)/f(r < 
500 kpc) computed as the integral of the surface brightness in 
the 0.5 - 5 keV band (observer frame). This quantity ranges be- 
tween and 1 and it represents the relative weight of the cen- 
tral surface brightness. Higher values of C are expected if a cool 
core is present. If the decrease in Zf e with redshift is associated 
to a decrease in the number of cool-core clusters for higher z, 
we would expect to observe a positive correlation between Zp e 
and C and a negative correlation between C and redshift. In 
Fig. [15] we plot Zf e as a function of C for our sample. We find 
that in our sample there is no correlation between metallicity 
and C with a Spearman's coefficient of r s = 0.02 (significance 
of ~ 0.2cr) nor one between C and redshift (r s = -0.1 1, a level 
of confidence of 0.8cr). 

The absence of strong correlations between C and iron 
abundance or between C and redshift suggests that the mix of 
cool cores and non cool cores over the redshift range studied 
in the present work cannot justify the observed evolution in the 
iron abundance. We caution, however, that a possible evolution 
of the occurrence of cool-core clusters at high redshift may still 
partially contribute to the observed evolution of Zf e . In other 
words, whether the observed evolution of Zp e is contributed 
entirely by the evolution of the mass of iron or is partially due 
to a redistribution of iron in the central regions of clusters is an 
open issue to be addressed with a proper and careful investiga- 
tion of the surface brightness of the high-z sample. 

A final check is provided by the scatter plot of Zf e versus 
Rext/Rvir, shown in Fig. [16] We do not detect any dependence 
of Zf e on the extraction radius adopted for the spectral analysis. 
In particular, we find that clusters with smaller extraction radii 
do not show higher Zf e values. 

Together with the previous discussion of possible selection 
biases of our sample, all these tests concur to indicate that the 
observed evolution of the iron abundance is a genuine signature 
of some physical processes associated with the production and 
release of iron into the ICM. 
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Fig. 15. Iron abundance plotted versus C = f(r < 
50 kpc)/ f(r < 500 kpc). Clusters within different redshift bins 
are coded with different symbols. 



This finding may be directly i nterpreted in terms of the cos- 
mic star formation history (see lEttorill2005) for high assume d 
values of the time delay of SNe la (see iDahlen et al.1 120041) . 
which are expected to be the ma in contributors of iron. Other 
works (e.g. lMannucci et al . 2006) argue that the data on (i) the 
evolution of the SN la rate with redshift, (ii) the enhancement 
of the SN la rate in radio-loud early type galaxies, and (iii) the 
dependence of the SN la rate on the colours of the parent galax- 
ies suggest the existence of two populations of progenitors for 
SN la. One population is expected to explode soon after the 
stellar birth on a time scale of 10 s years and can significantly 
pollute the ICM with iron at high redshift. The second popu- 
lation contributes to late enrichment with an exponential de- 
cay time of ~ 3 Gyrs. Following the argume nt described in 



Ettori (2005), we note that the rates adopted by lMannucci et al 
(2005) predict a flatter distribution of the iron abundance as 



function of redshift than the rates tabulated in Dahlen et al 



(2004), with a milder negative evolution still partially consis- 



I. Balestra et al.: Tracing the evolution in the iron content of the ICM 



15 



0.6 



0.2 



• 0.3<z<0.4 
■ 0.4<z<0.5 
A 0.5<z<0.65 
□ 0.65<z<1.0 
O 1.0<z<1.3 - 



O •* " A. a 



A o 
O °[ 

O A /\ 



•o ADi 




0.15 0.2 0.25 

R ext/ R vlr 



Fig. 16. Iron abundance plotted versus the ratio R ex t/Rvir- 
Shaded areas show the rms dispersion around the average iron 
abundance in four bins. 

tent with our measurements. Using detailed chemical evolution 
models, lLoewenstein (2006) has recently interpreted the sig- 
nificant enrichment of the ICM at z — 1 as direct evidence 
of prompt star formation in spheroids with a top-heavy IMF. 
Moreover, these models generally predict a significant increase 
in the iron abundance between z — and z - 1, in qualitative 
agreement with our results. 

Alternatively, the evolution in Z Fe occurring in the 5 Gyr 
spanned by our sample may be ascribed to some dynamical 
processes that transfer the metal-enriched gas from the inter- 
galactic medium of the cluster galaxies to t he hot phase of 



the ICM. Mechan isms such as ram pressure (Domainko et al 
2004 ICorall2006h or tidal stripping (e.g. Mur ante et alj|2004 



are currently being invest i gated with numerical simulations 



(see also Valdarnini 2003; Tornatore et al. 2004). The same 



mechanisms are often invoked to explain the morphological 
evolution of the cluster galaxies' population. 

5. Conclusions 

We have presented the spectral analysis of 56 clusters of galax- 
ies at intermediate-to-high redshifts observed by Chandra and 
XMM-Newton. This work improves our first analysis aimed at 
tracing the evolution of the iron content of the ICM out to z ^ 1 
(Paper I), by substantially extending the sample. The main re- 
sults of our work can be summarized as follows: 

- We determine the average ICM iron abundance with a 
~ 20% uncertainty at z > 1 (Z Fe = 0.27 ± 0.05 Z Q ), thus 
confirming the presence of a significant amount of iron in 
high-z clusters. Z Fe is constant above z — 0.5, the largest 
variations being measured at lower redshifts. 

- We find a significantly higher average iron abundance in 
clusters with kT < 5 keV, in agreement with trends mea- 
sured in local samples. For kT > 3 keV, Zp e scales with 
temperature as Z Fe (T) » 0.88 T^ )A1 . 



- We find significant evidence of a decrease in Zf e as a func- 
tion of redshift, which can be parametrized by a power law 
(Z Fe ) Z Fe (0)(l + z)~ a \ with Z Fe (0) =! 0.54 ± 0.04 and 
a z a 1.25 + 0.15. This implies an evolution of more than a 
factor of 2 from z - 0.4 to z = 1 .3. 

We carefully checked that the extrapolation towards z — 0.2 
of the measured trend, pointing to Z Fe - 0.5 Z , is consistent 
with the values measured within a radius r = 0.15 R V i r in lo- 
cal samples including a mix of cool-core and non cool-core 
clusters. We also investigated whether the observed evolution 
is driven by a negative evolution in the occurrence of cool-core 
clusters with strong metallicity gradients towards the center, 
but we do not find any clear evidence of this effect. We note, 
however, that a proper investigation of the thermal and chemi- 
cal properties of the central regions of high-z clusters is neces- 
sary to confirm whether the observed evolution by a factor of 
~ 2 between z = 0.4 and z — 1 .3 is due entirely to physical pro- 
cesses associated with the production and release of iron into 
the ICM, or partially associated with a redistribution of metals 
connected to the evolution of cool cores. 

Precise measurements of the metal content of clusters over 
large look-back times provide a useful fossil record for the 
past star formation history of cluster baryons. A significant iron 
abundance in the ICM up to z — 1.2 is consistent with a peak 
in star formation for proto-cluster regions occurring at redshift 
z — 4-5. On the other hand, a positive evolution of Z Fe with 
cosmic time in the last 5 Gyrs is expected on the basis of the 
observed cosmic star formation rate for a set of chemical en- 
richment models. Our data provide further constraints on the 
chemical evolution of cosmic baryons in the hot diffuse and 
cold phases. 
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Appendix A: Analysis of simulated X-ray 
spectra 

In this section we investigate the presence of a possible bias 
in the measure of Z/? c in our analysis procedure. In particular, 
we check whether unresolved gradients in the temperature or in 
the iron abundance distribution can affect the observed trends, 
paying particular attention to the low S/N regime of our spec- 
tra. We perform several simulations of spectra with different as- 
sumptions, as described in detail in the following subsections, 
and explore the possible conditions that can potentially affect 
the distribution of best-fit values of kT and, most important, of 
Zf e . The median of the distribution of best-fit values and the 
16% and 84% percentiles (corresponding to the lcr confidence 
level) will finally be compared to the input values of tempera- 
ture and metallicity. 

A. 1. Fitting bias in isothermal, constant metallicity, 
low S/N spectra 

A first simple test is to check the accuracy that we can achieve 
in recovering the input parameters of temperature and metallic- 
ity from spectra simulated with the typical S/N, temperatures 
and redshifts of clusters in our sample, under the assumption 
of a single-temperature mekal model. This may seem a redun- 
dant exercise. However a potential problem rises from the fact 
that upper limits on temperature are typically less constrained 
than lower limits. This effect increases at high temperatures, 
when the exponential cut off of the thermal spectrum shifts to 
energies for which the effective area of the detectors is low. As 
a consequence, if temperature best-fit values tend to be scat- 
tered upwards for low S/N, the estimated continuum may be 
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higher than the actual one, and therefore the measured equiv- 
alent width of the iron lines may be underestimated. On the 
other hand, variations in the best-fit temperatures affect the ion 
abundances too. Spectral simulations can be used to investigate 
how these aspects affect the measure of Zp e at low S/N as a 
function of redshift. 

In principle, the parameter space to explore is fairly wide: 
input metallicity, input temperature, redshifts, and S/N (mea- 
sured as S IN = C s l VGo/ + Cbcki where C s is the number of 
net counts from the source, C m the number of source plus back- 
ground counts, and Cj, c k the number of estimated background 
counts in the source area). For simplicity, we have decided to 
restrict the simulations to a few relevant cases. We chose an 
input metallicity of Zp e = 0.3 Z , redshift z = 0.4, and tem- 
peratures of kT = 3.5 and 7 keV. We took the observation of 
V 1416 as a template: the exposure time is 30 ks and the ex- 
traction radius 74". The simulations were performed for four 
different values of S/N (corresponding to 1850, 1000, 500, and 
200 net counts). We also run two simulations for z = 1. Each 
combination of parameters was simulated 1000 times. The sim- 
ulations were performed with XSPEC using a mekal model. 
We analyzed each simulated spectrum by adopting the same 
input model. In Table I A. 1 1 we list the median of the distribu- 
tions of the best- fit values for temperature and Zp e . Errors on 
the median are the 16% and 84% percentiles. 

The results are summarized in Fig. IA.ll We notice that, as 
expected, the low S/N spectra tend to have higher median best- 
fit temperatures compared to the input values. This translates 
into a slightly higher median of the best-fit Zf e than the input 
value, which is always Zp e = 0.3 Z for this set of simulations. 
However, for S /N < 20, the distribution of best-fit values is 
largely scattered around the input values. At higher redshifts, 
the situation for a given S/N improves slightly, since the expo- 
nential cutoff moves towards the most sensitive energy range 
of Chandra, and both the temperature and the abundance esti- 
mates are closer to the input values. We conclude that, under the 
assumption of a single temperature, single metallicity thermal 
plasma, the best-fit values of Zp e are not significantly biased, 
within the typical S/N and redshift range of our sample. 

A. 2. Fitting bias in two-temperature, 
constant-metallicity spectra 

Here we intend to investigate whether the presence of sub- 
stantial temperature structure can affect the measure of the 
iron abundance when the spectra are analyzed with a single- 
temperature mekal model, as adopted in this work. In some 
cases in our spectra, we are able to detect the presence of a 
temperature decrease towards the center; however, we are not 
able to perform a spatially-resolved spectral analysis for the 
large majority of clusters in our sample. 

Since there are no canonical, physically-motivated, multi- 
phase models of the ICM, we simply assume a double mekal 
model. Again, the free parameter space is wide, therefore here 
we only explore a few cases. In particular, we performed a set 
of simulations of a cluster with two temperature components 
(2 and 7 keV), with an emission measure ratio (EM, the nor- 
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Fig. A.l. Upper Panel: Median of the best-fit temperature dis- 
tribution as a function of S/N. Red solid circles refer to an 
input temperature value of 3.5 keV (dashed horizontal line), 
while blue solid squares to 7 keV (solid horizontal line). Empty 
squares are for z = 1. Lower and upper error bars correspond 
to the 16% and the 84% percentiles, respectively. Lower Panel: 
Median of the best-fit Zp e distribution as a function of S/N. 
The input value is always Zp e = 0.3 Z (horizontal dashed 
line). Symbols and error bars are as in the upper panel; ma- 
genta empty squares are for z — 1 . 



malization of the mekal model) of the cold to the hot com- 
ponent ranging from 0.3 to 0.75. The choice of the temper- 
ature range here agrees with the observational evidence that 
the minimum temperature in cool cores is always equal to or 
higher than one third of the ma ximum value of the hot com- 
ponent (see lPeterson et all2003l) . The iron abundance is set to 
0.3 Z in both components. We simulated 1000 spectra for each 
case at redshift z = 0.6 and z = 1. We took the observation of 
RX J0542 as a template: the exposure time is 50 ks and the ex- 
traction radius is 79". The results are listed in Table lA!2l and 
shown in Fig. IA.2I Obviously, for a given redshift, the best-fit 
temperature, which is always intermediate between the two in- 
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Table A.l. Results from the spectral simulations described in Appendix A.l. 
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Notes: a simulations identification number; b net number of counts; c signal-to-noise ratio; d input temperature; c input iron abundance; f 
redshift; g median values of the distribution of best-fit temperatures; h median values of the distribution of best-fit iron abundances. Lower 
and upper errors correspond to the 16% and 84% percentiles, respectively. 



put values, moves towards the low value for increasing values 
of the cold to hot EM ratio. But, for a given EM value, the me- 
dian of the best-fit temperature moves towards higher values 
for higher redshifts, since the cold components is redshifted 
out of the adopted energy range (0.6 - 8 keV). We note that, 
despite a two-temperature structure being fitted with a single- 
temperature, the fits are acceptable in the selected cases, due to 
the low number of total net counts (below 2000), representing 
the typical condition under which the presence of two temper- 
atures cannot be established from the spectral analysis. 

As for Zf e , the median of the distribution of the best-fit val- 
ues is slightly lower at z = 1 than at z = 0.6. The effect is a 
decrement of less than 30% up to z = 1 . Given the large disper- 
sion of the best-fit values, this effect is probably not playing a 
dominant role in our observed trend. However, a proper inves- 
tigation of the effect of a multi-temperature ICM must rely on a 
physical modelization of the ICM multiphase structure, which 
is presently missing. 

A. 3. Fitting bias in low S/N spectra with temperature 
and metallicity gradients 

We investigate the distribution of best-fit values for Zf e in the 
case of an ICM having significant temperature and metallic- 
ity structures at the same time. In particular, we investigate the 
most common case of higher Zf e associated with lower tem- 
perature components, as found in the cool cores. We repeated 
the simulations performed in Appendix A.2, by assigning an 
iron abundance of 0.6 Z to the cold component. The results 
are listed in Table lA.3l and shown in Fig. [A3] The best-fit tem- 
perature has the same behavior as described in Appendix A.2. 
For Zf c as well the situation is analogous to the previous case, 
with the median values slightly biased towards higher values. 
From a visual inspection of Figs. IA.2I and |A"31 we notice that 
the average best-fit values are lower for z — 1 compared with 
z - 0.6, by an amount on the order of 30%. On the basis of 
this result, we might expect that the presence of cool, metal- 
rich cores, could mimick the observed negative evolution from 
Z Fe a 0.4 Z at z ~ 0.3 to Z Fe a 0.2 Z at z ~ 1.3 with- 
out an effective decrease in the amount of iron in the central 



regions. However, this effect would explain only part of the ob- 
served trend even in the extreme case in which all the clusters 
in the sample had strong temperature and metallicity gradients. 
Therefore, we can conclude that the observed evolution of Zf e 
cannot be ascribed entirely to K-correction effects. We also no- 
tice that the limited effect of a central cold and metal-rich com- 
ponent in high-z clusters, also implies that the Zp e -temperature 
correlation cannot be simply explained by the occurrence of 
cool cores with temperature below 2 keV in clusters with virial 
temperatures kT < 5 keV. 

A.4. Fitting bias in isothermal spectra rich in 
a -element 

Finally, we checked whether a non-solar abundance ratio can 
affect the measure of Z Fe . In particular, we checked whether 
higher abundances of S, Si and O can artificially yield a higher 
Zp e . Therefore, we tried to recover the iron abundance with 
a mekal model, when the metallicity ratio among elements 
is higher than solar. The simulated spectra have the follow- 
ing input: Zf e = 0.3 and Z a = 1 or Z a - 2. We recall that 
these are extreme cases with respect to the abu ndances of a- 



elem ents observed in local X-ray clusters (see iTamura et al 
I2004T) . We simulated 1000 spectra of a z = 0.4 cluster with 
kT = 3.5 keV (2060 net counts expected) and with kT = 2 keV 
(2300 net counts expected), maximizing the presence of lines 
from a-elements in our energy range (lowest redshift and low- 
est temperature in our sample). We find that the temperatures 
are slightly higher that the input values in both cases (kT = 
3.8 ± 0.29 and 2.45 + 0.12 keV), while iron abundances are 
consistent (0.30 ± O.12Z and Z = 0.25^ Z for kT = 3.5 
and 2 keV, respectively). Therefore no bias is found in the mea- 
sure of Zf e , confirming the expectation that the iron abundance 
is mostly determined by the K-shell complex at 6.7 - 6.9 keV 
and is not affected by the L-shell complex below 2 keV, where 
the iron emission lines are blended to that of Si, O, and Mg. 
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Table A.2. Results from the spectral simulations of two temperature components, analyzed with a single-temperature mekal 
model, as described in Appendix A. 2. 
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Notes: a simulations identification number; b net number of counts of the hot and cold component; c emission measure ratio; d signal-to-noise 
ratio; e input temperature of the hot and cold component; f input iron abundance; ° redshift; h median values of the distribution of best-fit 
temperatures; 1 median values of the distribution of best-fit iron abundances. Lower and upper errors correspond to the 16% and 84% 
percentiles, respectively. 



Table A.3. Results from the spectral simulations of two components with different temperatures and iron abundances, analyzed 
with a single-temperature mekal model, as described in Appendix A.3. 
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Notes: a simulations identification number; b net number of counts of the hot and cold component; c emission measure ratio; d signal-to-noise 
ratio; e input temperature of the hot and cold component; f input iron abundance of the hot and cold component; g redshift; h median 
values of the distribution of best-fit temperatures; 1 median values of the distribution of best-fit iron abundances. Lower and upper errors 
correspond to the 16% and 84% percentiles, respectively. 
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Fig. A.3. Median of the best-fit Zf e distribution as a function 
of the ratio of the Cold and Hot Emission Measure, for simu- 
lations where the cold component has an iron abundance twice 
higher than the hot component (0.6 Z vs 0.3 Z ). Filled circles 
are for z — 0.6, while empty circles are for z = 1 . Lower and up- 
per error bars correspond to the 16% and the 84% percentiles, 
respectively. 



Fig. A.2. Upper Panel: median values of the best-fit tempera- 
ture distribution as a function of the ratio of the cold and hot 
emission measure. Input temperatures are 2 and 7 keV (hori- 
zontal dashed lines). Filled circles are for z = 0.6, while empty 
circles are for z = 1. Lower and upper error bars correspond to 
the 16% and 84% percentiles, respectively. Lower Panel: me- 
dian values of the best-fit Zp e distribution as a function of the 
ratio of the cold and hot emission measure. Symbols and errors 
as in the upper panel. 
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